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Abstract. Over the last few years enormous progress has been made in the numerical 
description of the inspiral and merger of binary black holes. A particular effort has gone 
into the modelling of the physical properties of the final black hole, namely its spin and recoil 
velocity, as these quantities have direct impact in astrophysics, cosmology and, of course, 
general relativity. As numerical-relativity calculations still remain computationally very 
expensive and cannot be used to investigate the complete space of possible parameters, semi- 
analytic approaches have been developed and shown to reproduce with very high precision the 
numerical results. I here collect and review these efforts, pointing out the relative strengths 
and weaknesses, and discuss which directions are more promising to further improve them. 
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1. Introduction 

Despite the almost unnatural simplicity with which the problem can be formulated (black 
holes are after all the simplest macroscopical objects we know), the final evolution of a binary 
system of black holes is an impressively complex problem to solve. At the same time, this very 
simple process plays a fundamental role in astrophysics, in cosmology, in gravitational-wave 
astronomy, and of course in general relativity. Recent progress in numerical relativity initiated 
by the works in refs. H]|2l[3l, have made it now possible to compute the different stages of 
the evolution, starting from the inspiral at large separations, for which post-Newtonian (PN) 
calculations provide an accurate description, through the highly relativistic merger, and finally 
to the ringdown. 

As long as the two black holes are not extremal and have masses which are not too different 
from each other (see however [|4) for simulations with a mass-ratio of 1 : 10), no major 
technical obstacle now prevents the solution of this problem in full generality and with an 
overall error which can be brought down to less than 1% (see ref. J3] |6] [7] [8] |4] for some 
recent examples). Yet, obtaining such a solution still requires a formidable computational 
power sustained over several days. Even for the simplest set of initial data, namely those 
considering two black holes in quasi-circular orbits, the space of parameters is too vast to 
be explored entirely through numerical-relativity calculations. Furthermore, many studies 
of astrophysical interest, such as many-body simulations of galaxy mergers, or hierarchical 
models of black-hole formation, span a statistically large space of parameters and are only 
remotely interested in the evolution of the system during the last few tens of orbits and much 
interested in determining the properties of the final black hole when the system is still widely 
separated. 

In order to accommodate these two distinct and contrasting needs, namely that of sampling the 
largest possible space of parameters and that of reducing the computational costs, a number of 
analytical or semi-analytic approaches have been developed over the last couple of years. In 
most of these approaches the inspiral and merger is considered as a process that takes, as input, 
two black holes of initial masses Mi, and spin vectors Si, S2 and produces, as output, a 
third black hole of mass Mfi n , spin Sf\ n and recoil velocity Vkick- Mathematically, therefore, 
one is searching for a mapping between the initial 7-dimensional space of parameters (i.e., the 
one containing the six spin components Si 2 and the mass ratio q = M2/M1) to two distinct 
3-dimensional spaces (i.e., the ones containing the three components of the final spin and of 
the recoil velocity Clearly this is a degenerate mapping (two different initial configurations 
can lead to the same final one) and it would seem a formidable task to accomplish given the 
highly nonlinear features of the few last orbits. Yet surprisingly, or perhaps not so surprisingly 
given the underlying simplicity of the problem, all of these studies have shown that the final 
spin vector and the final recoil velocity vector, can be estimated to remarkably good accuracy 
if the initial parameters are known f!Tl[l2l[T^[T4l[l5l[T6l[T7ll9l[T0l. 

This paper is dedicated to review the different approaches developed so far to tackle this 
problem, pointing out the relative strengths and weaknesses, and finally comparing them 
against the numerical data. The manuscript is organized as follows: in Sect. [2] I review the 
different methods that have been suggested to model the final spin, distinguishing those that 
are purely analytic from those that employ, with different amounts, the results of numerical 
simulations. Sect. [3] is instead dedicated to the specific approach developed at the AEI to 
obtain a simple and robust algebraic expression for the final spin vector. In Sect. [5] I also 

I The space of final parameters could in principle be larger if the final mass of the black hole were to be modelled. 
Two attempts to do so have been recently proposed (9l llOI (see also the discussion in Sec. [2), but this quantity is not 
usually modelled. 
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exploit the advantages of such analytic formula to explore those configurations that may be 
astrophysically more interesting and then validate the accuracy of the formula against all of the 
available data and across alternative approaches. Finally, Sect.|4]is dedicated to the modelling 
of the recoil velocity vector and a summary is presented of the different contributions and of 
the debate still present on one them. Section [5] will present the conclusions and highlight the 
prospects of future work in the modelling of the final state. 

As a final but important remark I note that all of the considerations made here apply to 
binary systems that inspiral from very large separations and hence through quasi-circular 
orbits. Such configurations are the ones more likely to occur astrophysically since any residual 
eccentricity is lost quickly by the gravitational -radiation reaction lfl8l . Furthermore, at least 
for nonspinning equal-mass black holes, recent work has shown that the final spin does not 
depend on the value of the eccentricity as long as the latter is not too large lfl9l . 

2. Modelling the final spin vector 

A number of analytical approaches have been developed over the years to determine the final 
spin from a binary black hole coalescence j20l [2T1 l22l l23l l24l . A particularly active line of 
research is the one that has exploited the motion of test particles in black hole spacetimes. 
A first attempt in this direction was made by Hughes and Blandford |25l . who assumed 
that the energy and angular momentum radiated in the merger and ringdown phase is much 
smaller than the one radiated during the inspiral. In this way they were able to express the 
conservation of energy and angular momentum when the two black holes reach the innermost 
stable circular orbit (ISCO) and hence compute the mass and spin of the final black hol^jj 
Of course, the notion of the ISCO is well defined in the test-particle limit (i.e., for Mi — > 0) 
and the approach of Hughes and Blandford reproduces this limit by construction and provides 
reasonable results for a large set of parameters. However, it also leads to incorrect final spins 
(i.e., Sf\ n /Ml n > 1) for very rapidly rotating black holes and even at very small mass ratios. 
Inspired in part by this work, Buonanno, Lehner, and Kidder [T2 j have more recently proposed 
a variant in which the angular momentum of the final black hole is assumed to be the sum of 
the individual spins and of the orbital angular momentum of a test particle at the ISCO of a 
Kerr black hole with the same spin parameter as that of the final black hole, i.e., 

M 2 a fin = MiM 2 L tp (|a fin |, fln ) + M x 2 &1 + A/f a 2 . (1) 

where ai t 2 =S\^IM\ 2 are the two dimensionless spin vectors (|ai, 2 | <G [0, 1]), 0& n is the 
angle of the final spin a^ n with respect to the direction of the initial orbital angular momentum 
of the binary, and L tp ( | dfi n | , #fin) is a vector whose direction is determined by assuming that 
the angle between L tp and the total spin S1+S2 maintains constant during the inspiral, and 
whose magnitude is the dimensionless angular momentum at the ISCO for a test particle 
moving around a black hole with spin |a,fi n | on an orbit with inclination 0g n with respect to 
the equatorial plane. 

Because L tp depends implicitly on the final spin afi n , equation ((TJ cannot be solved 
analytically and needs to be computed numerically. Furthermore, attention must be paid that 
L tp refers to a prograde (retrograde) orbit if the final spin is aligned (anti-aligned). 
Perhaps the most surprising aspect of this approach is that despite its simplicity, it is 
remarkably precise in a large portion of the space of parameters, with differences from the 
numerical-relativity results that can be as small as ~ 1% and of ~ 10% at most. As I will 
comment later on, I regard this result as an evidence that the dependence of the final spin on 

§ I note that although numerical simulations do not show evidence for the existence of an ISCO, the concept of an 
effective ISCO can nevertheless be useful for the construction of gravitational-wave templates [26 27 1. 
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the initial conditions is particularly simple and that, as a consequence, the mapping between 
the initial and final state can be accomplished with rather simple expressions. After all, L tp 
and ax, 2 have a limited functional excursion and are they linearly related in equation fl}. 
An important ingredient in the above recipe is "mass conservation", i.e., the assumption that 
the mass of the final black hole is the same as the sum of the two initial masses. Because 
gravitational radiation carries away part of the mass of the system, such an assumption is 
trivially incorrect, at least from a formal standpoint. However, although the assumption is 
false, the resulting approximation is not too severe. Numerical simulations reveal in fact that 
the mass radiated in gravitational waves is M Ta d/M = 1 — M Rn /M <5-7x 10~ 2 , with 
M = Mi + M2 being the total mass. Although this difference is of the order of the precision 
with which one wants to predict the final spin, a favourable series of cancellations reduces the 
impact of this otherwise unreasonable assumption. 

In the spirit of improving the accord with the simulations and in particular to produce more 
reliable predictions in the case of the merger of maximally spinning black holes, Kesden (9) 
has recently proposed an improvement over the approach in ref. Ifl2l which can be used for 
binaries with spins which are aligned or antialigned with the orbital angular momentum. Once 
again using the test-particle motion in a Kerr spacetime, Kesden has suggested that in addition 
to expression ([TJ being valid, the final mass Ma n for arbitrary mass ratios and initial spins can 
be computed as 

M Rn = M - fx[l - E(a Rn )} , (2) 

where /i = M\MijM is the reduced mass and E(a) is the energy per unit mass of a test 
particle on an equatorial orbit of a Kerr black hole with spin parameter a. This suggestion 
amounts to a minimal change in the calculation of the final spin, but provides a much more 
precise modelling of the final spin in the case of black holes with very large spins, i.e., with 
|<*i| - |a 2 | > 0.99. 

The resulting picture is particularly simple for equal-spin black holes, i.e., with \a\\ = | ci2 1 = 
a, (see Fig. 2 of ref. [9|). In this case the test-particle approach suggests that for a < 0.9916, 
the final spin afi n is a function increasing monotonically with the symmetric mass ratio 
v = M\Mil (M1+M2) 2 = q/(l+q) 2 , so that equal-mass mergers (i.e., with v = l/4)arethe 
most efficient ones in spinning up the final black hole. However, for values of a > 0.9916 the 
function aa n (v) is no longer monotonically increasing with v but has a maximum at smaller 
and smaller values of v, so that for these spins a suitable symmetric mass ratio v 7^ 1/4 
maximizes the final spin. This behaviour continues for increasingly large values of a and up 
to a ~ 0.9987, when a Rn (v = 1/4) = a; as aresult, for a > 0.9987 equal-mass mergers spin 
down black holes rather than spinning them up, i.e., a& n (v = 1/4) < a for a > 0.9987. 
Finally, when the initial black holes are maximally spinning, i.e., a = 1, the maximum 
of afi n (i/) is at the extreme mass-ratio limit (EMRL), i.e., at v = 0, thus making a^iy) 
monotonically decreasing with 4jJ As a direct consequence, maximally spinning black holes 
cannot be produced by the merger of two black holes, the only exception being offered by the 
EMRL. Rather, the merger of maximally spinning equal-mass black holes yields a black hole 
with o fin ~ 0.998cfiJ 

An alternative and distinct line of research for the calculation of the final spin has instead 
exploited numerical-relativity calculations as soon as the first successful evolutions were 

I Wald has shown in ref. 1281 that the cosmic censorship hypothesis holds also in the neighbourhood of the EMRL, 
i.e., that a^ly) < 1 for v ~ 0. 

K A different astrophysical route to the production of maximally spinning black holes is of course that of accretion 
and this was shown by Bardeen in ref. 1291 . The arguments made there apply also here if one is considering a binary 
coalescence in the neighborhood of the EMRL, i.e., for v ~ 0. 
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possible [30, 31 , 32 1 . While these first approaches essentially represented fits to the numerical 
data (which were initially scarce) they provided the first accurate description of the problem, 
albeit in a small region of the space of parameters. Subsequent work focussed instead on 
the derivation of analytic expressions which would model the numerical-relativity data but 
also exploit as much information as possible either from perturbative studies, or from the 
symmetries of the system when this is in the weak-field limit ifTTl IT3l IT4l IT31 IT61 ITtTI . In this 
sense, these approaches are not blind fits of the data, but, rather, use the numerical-relativity 
data to construct a physically consistent and mathematically accurate modelling of the final 
spin. 

The common ground shared by the approaches in this second class is given by the assumption 
that the final spin, when seen as the function afi n = ag n (ai, 02, v), can be expressed as a 
Taylor expansion around a\ = 0,% = v = 0. Given that |ai a| < 1, this may seem as a 
mathematically reasonable assumption and the expectation that the series is convergent over 
the whole space of parameters as a legitimate one. However, this remains an assumption and 
more work is needed to make this mathematical assumption also a physically reasonable one. 
Sharing this assumption, then different routes are chosen to constrain the coefficients and these 
may invoke more mathematically -based considerations, as those proposed in ref. lfT3lfT4l[T0l . 
or more physically -based considerations, as those proposed in refs. IfTTl [161 fTTj . Both 
approaches often reach the same conclusions, but there is little doubt that the systematic 
approach proposed in refs. lfT3l [141 offers the advantage of being more easily expandable 
to capture higher-order effects. 

Here, however, I will concentrate on reviewing the approach which, with a minimal 
number (4) of physically reasonable assumptions and with a minimal number (5) of free 
coefficients to be fixed from the numerical data, leads to a formula that can model generic 
initial spin configurations and mass ratios, thus covering all of the 7-dimensional space of 
parameters IfTTl [161 [171 . As I will show in Sect. [3] besides being simple and physically 
motivated, it is also remarkably accurate in reproducing the final spin of more than 150 
simulations. 

In essence, the approach developed in refs. IfTTl IT6l fTTl (hereafter simply the "AEI formula") 
amounts to considering the dimensionless spin vector of the final black hole as given by 
the sum of the two initial spins and of a "third" vector parallel to the initial orbital angular 
momentum when the binaries are widely separated. This "third" vector is an intrinsic 
"property" of the binary (it will be shown below that this is essentially the orbital angular 
momentum not radiated), thus depending on the initial spin vectors and on the black holes 
mass ratio, but not on the initial separation. The formula for the final spin then simply 
describes the properties of this vector (so far only of its length) in terms of the initial 
parameters of the binary and of a set of coefficients to be determined from a comparison 
with numerical simulations. 

Let us now consider in more detail how to derive such a formula. As mentioned in the 
Introduction, four assumptions are needed in order to make the problem tractable analytically 
and these are listed in what follows. 

(i) The mass radiated to gravitational waves M ra d can be neglected i.e., Ma n = Mi + Mi. 

As mentioned before, while this assumption is certainly not correct, its influence on the 
overall accuracy of the prediction is small, with the possible exception of binaries with very 
small mass ratios and with rapidly spinning black holes |9). Work is in progress to relax this 
assumption and will be presented in a forthcoming paper 11331 . 

( ii) At a sufficiently large but finite initial separation the final spin vector Sft n can be 
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approximated as the sum of the two initial spin vectors and of a third vector £ 

s fi „ = s 1 + s 2 + J, 

When viewed as expressing the conservation of the total angular momentum, eq. also 
defines the vector £ as the difference between the orbital angular momentum when the 
binary is widely separated L, and the angular momentum radiated until the merger J ra d, 
i.e., £=L—J la< i. Stated differently, the vector £ measures the orbital angular momentum that 
cannot be radiated and, assuming the "effective" ISCO as the radial separation at which the 
system essentially stops radiating, £ can then be thought of as the angular momentum of the 
binary at such ISCO. 

( Hi) The vector £ is parallel to L. 

This assumption (which is made also in ref. |[T2"1 ) essentially enforces that the component 
of the final spin Sf\ n in the orbital plane equals the one of the total initial spin S1+S2 in 
that plane. This is correct when Si = —S2 and q = 1 [this can be seen from the PN 
equations at 2.5 order], or by equatorial symmetry when the spins are aligned with L, or 
when Si =S2 = (also these cases can be seen from the PN equations). However, for 
more general configurations one expects that £ will also have a component orthogonal to L 
as a result, for instance, of spin-orbit or spin-spin couplings, which will produce in general 
a precession of £ . In practice, the component of £ orthogonal to L will correspond to the 
angular momentum J^ d radiated in a plane orthogonal to L, with a resulting error in the 
estimate of | £ | which is- | J^ d \ 2 /\£ | 2 ~ | J^ d \ 2 /(2y/3M 1 M 2 ) 2 . A direct consequence of 
assumption ( Hi) is a systematic overestimate of the final spin component along the direction 
of the orbital angular momentum a{| n and, consequently, a systematic overestimate of < 2% 
in the cosine of the inclination angle 6*fi n . In ref. iflOl this assumption was relaxed and the 
modelling of a\ n improved; however, a comparison with the available data will show that 
the error introduced by this assumption is not important for the estimate of the final spin 
magnitude |a|fi n (see Fig.[3]and the discussion in Sect. [4]). 

(iv) When the initial spin vectors are equal and opposite (S± = —S2) and the masses are 
equal (q = 1), the spin of the final black hole is the same as for nonspinning binaries. 

Stated differently, equal-mass binaries with equal and opposite-spins behave as nonspinning 
binaries, at least as far as the properties of the final black hole are concerned. This 
condition is met by the leading-order contributions to the spin-orbit and spin-spin point- 
particle Hamiltonians and by the spin-induced radiation flux l34l [22], but it is difficult to 
justify in full general relativity. Here it is introduced essentially to break a possible degeneracy 
in the formula (see below) and it reflects the expectation that if the spins are the same and 
opposite (independently on whether they are aligned/antialigned with the orbital angular 
momentum), their contributions to the final spin cancel (at least to the precision of interest 
here) for equal-mass binaries. Besides being physically reasonable, this expectation is met by 
all of the simulations performed to date, both for spins aligned with L lfTTl[T6l and orthogonal 
to L ll35l [TBI [lOl . A similar assumption is also made, although not explicitly, in ref. Ifl2l . 

Using these assumptions it is now possible to derive the analytic expression for the final spin. 
Let us start by expressing the vector relation (O as 

atin = ^ (ai + a 2 q 2 + iq) , (3) 

where a^ n =Sa n /M 2 [cf. assumption f/J], £= £/ (M1M2), ai,2 = S\.2/M 2 2 , and its norm 
is then given by 
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a fin 



(1+9) 2 



|ai| 2 + |a 2 |V + 2|a 2 ||ai| 9 2 cosa + 



1/2 



2(|a 1 |cos/3+|a 2 | 9 2 cos7)|^|g+|£|VJ ' , (4) 

where the three vector cosine a, (3 and 7 are defined by 

cosa = ai-a 2 , cos/3 = ai-£, cos7 = a 2 -^. (5) 

Because ai l2 || Si )2 and £ \\ L, the angles a, /? and 7 are also those between the initial spin 
vectors and the initial orbital angular momentum. Indeed, it is to exploit this simplification 
that assumption (Hi) was introduced and thus to replace &i >2 with Si. 2 and £ with L in (0. 
Note that in this line of arguments I am making the underlying assumption that a, (3 and 7 are 
well-defined and hence that the initial separation of the two black holes is sufficiently large 
that the the changes in these angles due to precession is small. The error introduced in this 
way in the measure of cos a, cos (3 and cos 7 are expected to be of the order of ~ |^ r 4dl/l^l- 
All that is needed at this point is to "measure" \£\ by exploiting the results of numerical- 
relativity simulations and, in particular, those for binaries with spins parallel and aligned 
{i.e., a = (3 = 7 = 0), for spins parallel and antialigned spins (i.e., a = 0, (3 = 7 = n), and 
for antiparallel spins which are aligned or antialigned (a = f3 = tt, 7 = or a = 7 = ir, 
= 0). 

As mentioned before, the matching with the numerical data is not unique but this degeneracy 
can be broken if assumption (iv) is adopted. Doing this and requiring in addition that 
\£\ depends linearly on cos a, cos/3 and cos 7, one finally obtains that the modulo of the 
dimensionless, non-radiated angular momentum is given by (see ref. JT7] for a more detailed 
derivation of this expression) 

\ £ \ = TJ^-m a i| 2 + MV + 2|ai||a 2 |g 2 C osa) + 

s 5 v + to + 2 \ (| ai | cos/3+ | a2 | g 2 ) + 2V3 + t 2 is + t 3 is 2 . (6) 

1 + <r / 

The 5 coefficients S4, S5, to, t 2 , £3 appearing in (O were determined by fitting a large number 
of binaries having equal-mass and unequal but aligned/antialigned spins ifTTl or unequal-mass 
but equal aligned/antialigned spins lfl6l . As a result of this fitting process, the coefficients 
have been determined to be 

s 4 ^-0.129, s 5 ~ -0.384, 

i ~ -2.686, t 2 ~ -3.454, 

t 3 ~ 2.353. (7) 

Expressions and © fully and accurately determine the modulo of the final spin under 
generic spin configurations and mass-ratio, but they do not yet provide information about the 
different components of the final spin vector. Such information can be extracted from the 
angle between the final spin vector and the initial orbital angular momentum #fj n , which can 
be easily calculated from |agn|. Because of assumption (Hi), the component of the final spin 
in the direction of L is [cf. eq. (O] 

H ? |ai|cosff+ |a 2 |g 2 cos 7+ \£\q 

Kn = a fin • t = {l+q) 2 ' ( } 

so that cos0fi n = aj[ n /|afm|, and the component orthogonal to the initial orbital angular 
momentum is a^ n = ag n sin0g n . Note that by construction our approach does not further 
decompose the perpendicular component of the spin a^ n . 
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Before moving on to a more detailed exploration in Sect. |4] of the physical implications of 
expressions (0]i and (O, it is useful to remark that the approach reviewed here is intrinsically 
approximate. Although the comparison between the predictions of the formula and all of the 
available numerical data has revealed differences of a few percent at most, attention should be 
paid to the fact that the coefficients have been constrained with a rather small and statistically 
non-homogeneous set of configurations. Most importantly, however, the formula discussed 
above can be improved in a number of different ways which I list below: (a) by reducing the 
X 2 of the fitting coefficients as new simulations are carried out; (b) by using fitting functions 
that are of higher-order than those used so far [cf. expressions (TTOb and (0] in the spirit of 
the systematic spin expansion suggested in ref. lfT3l IT41 ; (c) by estimating J^ d through PN 
expressions or by measuring it via numerical simulations; (d) by introducing an independent 
modelling of the mass radiated in gravitational waves M ra( j. Work is in progress to include 
some of these possible improvements and it will be presented in l33ll . 

3. Exploring the AEI spin formula 

In what follows I discuss in some detail the predictions of the AEI spin formula for some 
simple cases and highlight how to extract interesting physical considerations. 

3.1. Unequal mass, Aligned/Antialigned Equal spins 

If the black holes have unequal mass but spins that are equal, parallel and aligned/antialigned 
with the orbital angular momentum, i.e., \ai\ = |oa| = a, a = ; (3 = 7 = 0, tt, the 
prediction for the final spin is given by the simple expression lfl6l 

afi n (a, v) = a cos /3 (l + s^av + s 5 is 2 + t v) + 2v3i/ + t 2 v 2 + t 3 v 3 . (9) 

where cos/3 = ±1 for aligned/antialigned spins. Note that since the coefficients in (0 
are determined by fits to the numerical data and the latter is scarcely represented at very 
large spins, the predictions of expression (0 for nearly maximal black holes are essentially 
extrapolations and are therefore accurate to a few percent at most. As an example, when 
a = 1, the formula (0 is a non-monotonic function with maximum a.fi n ~ 1.029 for 
v ~ 0.093; this clearly is an artifact of the extrapolation. 

The global behaviour of the final spin for unequal-mass and aligned/antialigned equal-spin 
binaries is summarized in Fig. [T] which shows the functional dependence of expression (0 
on the symmetric mass ratio and on the initial spins. Squares refer to numerical-relativity 
values as reported in refs. ||36l|33|38][l]3[l3J[39]|4), while circles to the EMRL constraints. 
A number of interesting considerations can now be made: 

(a) Using (0 it is possible to estimate that the minimum and maximum final spins for an 
equal-mass binary are afi n = 0.3502 ±0.03 and afi n = 0.9590 ±0.03, respectively. While the 
value for the maximum spin is most likely underestimated (indeed Kesden predicts for this 
configuration aa n = 0.9987 |9|) the minimum value is expected to be much more accurate 
than the estimate in ref. [0, which tends to underestimate the final spin for a < —0.3 (cf. Fig. 
5 of ref. 0. 

(b) Using (0 it is straightforward to determine the conditions under which the merger will lead 
to a final Schwarzschild black hole. In practice this amounts to requiring aa n (a, v) — and 
this curve is shown as a (red) dashed line in Fig.Q] Several numerical simulations have been 
carried out to validate this condition lfl6l [39l and all of them have shown to produce black 
holes with afi n < 0.01 (cf. squares in Fig.Q]with v ~ 0.16). Overall, the behaviour captured 
by expression (0 shows that in order to produce a nonspinning black hole it is necessary to 
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numerical data □ 




Figure 1. Global dependence of the final spin on the symmetric mass ratio and on the initial 
spins as predicted by expression {9) for equal-mass, aligned/antialigned equal-spin binaries. 
Squares refer to numerical-relativity values as reported in refs. [36 37 38 16 15 1391 l4l. 
while circles to the EMRL constraints. Indicated with a (red) dashed line is the locus of points 
leading to a Schwarzschild black hole (i.e., ag n = 0), while (green) solid lines mark the region 
leading to a "spin-flip" (i.e., aa n a < 0). 

have unequal-masses (the largest possible mass ratio is v ~ 0.18) and spins antialigned with 
the orbital angular momentum to cancel the contribution of the orbital angular momentum to 
the total one. 

(c) Using © it is also straightforward to determine the conditions under which the merger 
will lead to a "spin-flip ", namely when the newly formed black hole will spin in the direction 
opposite to that of the two initial black holes. Mathematically this is equivalent to determine 
the region in the plane (a, v) such that afj n (a, v) a < and it is shown in Fig.Q]as limited by 
the (green) solid lines and by the (red) dashed line. Overall, it is clear that a spin-flip can take 
place only for very large mass ratios if the black holes are initially rapidly spinning and that 
small mass ratios will lead to a spin-flip only for binaries with very small spins. 

(d) Finally, using (O it is also possible to determine the conditions under which the merger 
will lead to a final black hole with the same spin as the initial ones. This amounts to requiring 
that afj n (a, v) — a = and only a very small portion of the (a, v) plane does satisfy this 
condition (cf. Fig. 5 of ref. lfl6ll ). For equal-mass binaries, for instance, the critical value is 
Qcrit ^ 0.946 and no spin-down is possible for v < 0.192. Because of the minuteness of the 
region for which a^ n < a, black holes from aligned-spins binaries are typically spun-up by 
mergers. As it will be shown also in the following Section, this statement is true also for other 
configurations and is probably true in general. 

3.2. Equal-mass, Aligned/ 'Antialigned Unequal spins 

Equally interesting is to consider the prediction for the final spin in the case in which the 
initial black holes have equal mass but unequal spins that are either parallel or antiparallel 
to the orbital angular momentum, i.e., for q = 1 and a — 0, ir ; (3 — 0, tt ; 7 = 0, it. Setting 
2\ai I cos j3 = ai + a,2 in expression (0]i we obtain the simple expression for the final spin in 
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Ofin (01,02) =Pa+Pi{ai + a 2 )+p 2 (a 1 + a 2 ) 2 , (10) 
where the coefficients po,Pi and P2 are given by 

P2 = £± ~ -0.0081 . (11) 
16 

Note that the coefficients po,Pi,P2 and S4, S5, to, t 2 , £3 were obtained through independent 
fits of two distinct data sets. The fact they satisfy the conditions (fTTT i within the expected 
error-bars is an important consistency check. 

When seen as a power series of the initial spins, expression ( TTOb suggests an interesting 
physical interpretation. Its zeroth-order term, po, can be associated with the (dimensionless) 
orbital angular momentum not radiated in gravitational waves and thus amounting to ~ 70% 
of the final spin at most. Interestingly, the value for po is in very good agreement with what is 
possibly the most accurate measurement of the final spin from this configuration and that has 
been estimated to be a fin = 0.68646 ± 0.00004 ||40). Similarly, the first-order term in ( [Tol l, 
Px, can be seen as the contributions from the initial spins and from the spin-orbit coupling, 
amounting to ~ 30% of the final spin at most. Finally, the second-order term, p 2 , can be 
seen as accounting for the spin-spin coupling, with a contribution to the final spin which is of 
— 4% at most. 

Another interesting consideration is possible for equal-mass binaries having spins that are 
equal and antiparallel, i.e., q = 1, a\ = — a 2 . In this case, expressions (0]i and © reduce to 

i» 1 l £ l ^o. h x t3 

afin = — = ■ (1^) 

1 1 4 2 16 64 V ; 

Because for equal-mass black holes which are either nonspinning or have equal and opposite 

spins, the vector £ does not depend on the initial spins, expression (fT2l states that |£|AY| n /4 = 

|£|M 2 /4 = |£|MiM2 is, for such systems, the orbital angular momentum at the effective 

ISCO. We can take this a step further and conjecture that \£\MiAl2 is the series expansion 

of the dimensionless orbital angular momentum at the ISCO also for unequal-mass binaries 

which are either nonspinning or with equal and opposite spins. The zeroth-order term of this 

series (namely, the term 2\fiM\M 2 ) is exactly the one predicted from the EMRL. 



3.3. Generic (misaligned) binaries: unequal mass, unequal spins 

When the binaries are generic, namely when the initial spins are oriented in generic directions 
and the two masses are different, the spin formula (H}-© does not reduce to a simple 
expression and the analysis of the physical implications becomes more complex. 
Much more challenging is also the numerical solution in these cases, partly because they are 
computationally more expensive (no symmetries can be exploited to reduce the computational 
domain), and partly because the evolutions start at a finite separation which does not account 
for the earlier evolution of the orbital angular momentum vector and of the spins (both of 
which precess). In addition, because the final spin is oriented in directions which are in 
principle arbitrarily far from the main coordinate lines, the calculation of the inclination angle 
from the properties of the final apparent horizon is often non-trivial and suitable definitions 
need to be introduced (see, e.g., 1361 ). 

Notwithstanding this, accurate numerical simulations have been reported by several groups 
and Fig.|2]reports data for 75 of such binaries. Using a "dummy index" to order the different 
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dummy index 

Figure 2. Top panel: Comparison between the final spin values | ag n | obtained via numerical 
simulations and via the analytic formulas of AEI 1171 (red squares) and of FAU 1101 (green 
triangles). The numerical data has been published in ref. 1101 131 1 ETl 1 171 1421 171 (see main 
text for a more detailed description).Boffom panel: Residuals |afl n | — |a pr cd.l f° r tne AEI 
formulas (red squares) and for the FAU formula (green triangles). 



binaries, here is where the different results were presented: the data from 1 to 3 1 was presented 
in ref. flH, the data from 32 to 34 in ref. EHUD, the data from 35 to 37 in ref. H3, the data 
from 38 to 45 in ref. 11421 . and finally the data from 46 to 75 was presented in ref. Q. The 
top panel of Fig. [2] in particular, reports a comparison between the final spin values |afi n | 
obtained via numerical simulations and via the analytic formulas of AEI ifTTl (red squares) 
and of FAU iflOl (green triangles). As mentioned before, the latter has been proposed as a 
suitable application of the systematic expansion suggested in |13| when this is truncated at 
first order. 

Because derived from a generic expansion and not from some assumptions (as for the 
AEI formula), the the FAU formula offers the important improvement of going beyond the 
restrictions imposed by our assumption ( Hi) and, in particular, models the changes of the final 
spin orthogonal to the orbital angular momentum. I recall, in fact, that as a consequence of 
assumption ( Hi), the final spin component in the orbital plane, a^ n does not change in the 
AEI formula and it is simply that of the initial total spin Si+S^- Tichy and Marronetti IflOl . 
instead, can model this change both via a contraction (of about 70%) and a rotation (which is 
rather small, with a rotation angle which is less than 9 degrees). This leads to a better estimate 
of the inclination ang id3 but not necessarily of the final spin modulo. 

This is shown in the bottom panel of Fig. [2] where I report the residuals |<2fi n | — |<J P red.| for 
the AEI formulas (red squares) and for the FAU formula (green triangles). While the two 
approaches are overall comparable, the FAU one has larger residuals in general. This may 
be combined result of having used a small set of simulations to constrain the coefficients (the 
FAU formula has in fact used only the data from index 1 to index 31 and for these points 

+ To fix the ideas and using as a reference the data published in ref. (7J, the FAU formula has an average error on 
cos 6fj n of ~ 0.1%, while the AEI formula has an average error of ~ 1.7%. 
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Figure 3. Absolute residuals between the numerically computed value for the final spin and 
the one predicted by either the AEI formula (red squares) or the FAU formula (green triangles), 
i.e., | |<Xfi n | — |q pre < j | |. The numerical data is the one published in refs I36ll371f38l lll|[l6l 
[l5][^|3[T0][^Ef][i7]|42][7]|4O]|4l and is shown simply in terms of the initial spin moduli 
l tt |l,2- 

the residuals are clearly smaller) and of truncating the series expansion at first order (the AEI 
formula has exploited a third-order expansion in the spins and mass ratio). 
This conclusion applies also when comparing the two formulas against all of the available 
numerical data, that is all of the data reported for the 153 simulations presented in refs. l36l 
G3[3S[IT][T3[l3J[39]H This large bulk of data is shown in a 

very synthetic form in Fig. [3] where I report the absolute residuals between the numerically 
computed value for the final spin and the one predicted by either the AEI formula (red squares) 
or the FAU formula (green triangles), i.e., \ \aa n \ — |a pro( j.| |. Also indicated with constant 
planes are the arithmetic average of the residuals for the AEI formula (i.e., 0.775 x 10~ 2 ) and 
the arithmetic average of the residuals for the FAU formula (i.e., 1.011 x 10 -2 ). 
A few remarks are worth making. Firstly, the two formulas yield very similar results and 
while the AEI one is better on average, it also has the largest residual and all for the data 
from ref. l43l . Secondly, configurations with very high spins are those that produce the 
largest residuals and this is partly due to the difficulty of having accurate calculations in those 
regimes and, in part, to the fact that the coefficients in AEI's formula were determined from 
simulations with moderate spins. Finally, and most importantly, the data refers to simulations 
from different codes, different numerical setups and hence different truncation errors. It is 
therefore remarkable that the scatter is so small, with the largest majority of residuals being 
below 1%. Stated differently, Fig.[3]can be interpreted as an absolute variance of the precision 
of present numerical-relativity calculations. 

4. Modelling the final recoil vector 

As mentioned in the Introduction, the final spin is not the only quantity which is possible 
and useful to model from binary black-hole coalescences. Rather, the gravitational recoil 
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is a common feature of binary black-hole coalescences and with multiple consequences. 
I recall that a binary with unequal masses or unequal spins radiates gravitational energy 
asymmetrically and the linear momentum lost at any given time is not compensated by the loss 
after one period. As a result, during the inspiral the linear momentum vector of the asymmetric 
binary changes continuously with time increasing and spiraling outwards up until the two 
black holes merge and the final black hole conserves the linear momentum at the time of the 
merger. This recoil velocity, or simply "kick", was predicted well before numerical relativity 
simulations were able to prove its existence and measure its magnitude l44l [431 . Analytic 
estimates of the kick velocities have also been available for some time K6l [47l l48l [491 IBTfl . 
but because the largest part of the system's acceleration is generated in the final parts of the 
last orbit, when the motion is relativistic and the dynamics highly nonlinear, fully relativistic 
calculations are necessary in order to determine the kick accurately. 

Indeed, modelling the recoil velocity has been one of the most exciting results of recent 
numerical-relativity calculations, and these have had a deep impact not only in general 
relativity but also on astrophysics and cosmology. Discussing in detail such implications 
would be very interesting but impossible with the space restrictions of this contribution. 
Hereafter, therefore, I will only briefly review the main results obtained over the last couple of 
years and provide a snapshot of the present understanding of the recoil velocity from binary 
black-hole coalescences. 

Computations of the recoil velocity were made soon after the first successful binary black 
hole evolutions and the systems which were studied first were unequal-mass systems with 
moderate mass ratios [51 , 52]. These initial studies were quite limited in the mass ratios they 
were able to consider and were further extended in ref. l53l . where a large number of models 
with mass ratios between 0.25 and 1.0 was explored. A reasonably consistent picture has 
emerged from these simulations and it is one in which the recoil velocities generated by mass 
asymmetries are rather accurately bounded above by a maximum velocity of ~ 175km/s, 
which is attained for v ~ 0.195, or q ~ 0.36. Although very small ratios (i.e., q < 0.01) still 
represent the frontier of what is presently possible to compute with standard finite-difference 
codes, the recent work of ref. |4| with a mass ratio q = 0.1 suggests that no major surprise 
should be expected in the functional dependence of the recoil velocity for q < 0.1. 
As mentioned above, a recoil velocity is generated also when the asymmetry is in the spins 
rather than in the masses. A number of groups have therefore looked into the kick produced 
by equal-mass binaries having, however, unequal spins 14311341 1331 HTl 131 1361 . In this context, 
particular attention has been paid to binary black hole systems with equal-mass and spins 
aligned with the orbital one, since systems of this type may represent a preferred end-state of 
the binary evolution. PN studies in vacuum have shown, in fact, that the gravitational spin- 
orbit coupling has a tendency to align the spins when they are initially close to the orbital 
one 1371 . Furthermore, if the binary evolves in a disc, as expected for supermassive black 
holes, the matter can exert a torque tending to align the spins l58l . Overall, the results 
of several groups now agree that the recoil velocities generated from the merger of equal- 
mass, aligned/antialigned unequal-spins are rather accurately bounded above by a maximum 
velocity of ~ 450 km/s, which is attained for |ai| = — |ot2 1 = 1 (cf. Fig. 2 of ref. ifTTl ). Such 
velocities are large enough to eject the merged black hole from the center or a dwarf elliptic 
or spheroidal galaxy, and which have escape velocities < 300 km/s 1391 . 
Since the recoil is the result of an asymmetry in the binary system, it does not come as a 
surprise that the largest recoils are those produced in the most asymmetric systems, namely 
in those binaries where the spins are antiparallel and lie in the orbital plane. These are indeed 
the configurations leading to the "superkicks", i.e., extremely large kicks which have been 
computed to be of the order of 2000 — 2500 km/ s and then extrapolated to yield velocities as 
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large as 4000 km/s for the maximally spinning case BT1 l60l |6T1 [35l l62l . Such velocities are 
so large to be able to eject the merged black hole from the center of even massive elliptical 
galaxies (whose escape velocities are < 2000 km/s l59l ) and may represent a problem if this 
process does indeed take place frequently during the merger of massive galaxies. 
To clarify the picture of what contributes to what in the final recoil, it is useful to adopt the 
phenomenological splitting first proposed in ref. ATI , and which distinguishes the various 
contributions to the recoil velocity as coming from two different sources of asymmetry in 
the system, namely the one relative to mass-asymmetry v M and the one relative to the 
spin-asymmetry v s . The latter is then further decomposed in the component contained in 
the orbital plane (which is produced by the component of the spins parallel to the orbital 
angular momentum) and in the component parallel to the orbital angular momentum (which 
is produced by the component of the spins orthogonal to the orbital angular momentum). 
Selecting a Cartesian coordinate system in which the unit vector in the z-direction e% is 
aligned with the orbital angular momentum L, the above decomposition leads to a generic 
expression 

Vkick = v M e! + uj-(cos£ei + sin£e 2 ) + u|e 3 . (13) 

It should be noted that expression ( TT3l is essentially heuristic although the different 
contributions can have approximate expressions as derived from PN considerations. Over 
the last couple of years a considerable effort has been invested in modelling the different 
contributions which can now be expressed as follows. The mass-asymmetry contribution is 
given by J46]|53) 

v M =Av 2 (l + Bv)VT=4^, (14) 
the "in-plane", spin-asymmetry is given by l43l l54ll54ll4"T1 131 l56lfTT1 




while the "off -plane", spin-asymmetry is given by l60ll6Tl[35ll63ll64ll 

u " =S (TT^) l°i -«°a|cos(e-e ). (16) 

Here, A ~ 1.2 X 10 4 km/s, B ~ -0.93, C ~ 7040 km/s, D ~ 1460 km/s, and 
E ~ 6.0 x 10 4 km/s. Note that the angle is defined "a-posteriori" as the angle between the 
in-plane component of the vector A= (Mi + M2){S\/Mi— S2/M2) and the direction along 
which the two black holes approach each other at the merger, while Oo is a simple constant 
offset. 

A number of considerations are in order. Firstly, it should be noted that oc a| and 
viceversa 1A oc aj-, thus suggesting that the unbalanced gravitational-wave emission {i.e., the 
one responsible for the changes in of linear momentum vector) is always orthogonal to the 
black holes spins. Secondly, the in-plane spin-asymmetry component given by expression 
( TTol ) shows that the largest contributions come from binaries with anti-parallel initial spins, 
i.e., a\ = —a\ and that, in this case, the quadratic component is identically zero. Thirdly, 
the mass dependence of the off-plane spin-asymmetry component given by expression ( fTol l 
is presently under debate and indeed the work carried out in ref. l65l suggests instead a 
dependence of the type 

v s = E 7TT—\ (? a 2 cos (^ - $ 2) - ai cos(0i - $1)) . (17) 
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where $12 are constant offsets and </>i.2 are the angles between a^ 2 and some reference 
direction in the equatorial plane. The most conspicuous and important difference between 
expressions (TToT l and (fTTJi is in the dependence on the symmetric mass ratio (i.e., v vs 
and this is not a minor difference since it has important astrophysical consequences (see the 
discussion in ref. Il65l ). It is presently unclear which of the two expressions will find further 
confirmation and work is in progress ll66l to resolve this debate through an independent set of 
simulations and a careful modelling of the integration constant. The latter, I recall, accounts 
for the nonzero initial linear momentum the spacetime has at the time the simulation is started. 
An accurate determination of the integration constant has played an important role in the 
determination of the second-order correction in the contribution for vj: (see discussion in 
ref. 0). 

5. Conclusions 

The determination of the final spin and recoil velocity from the from the knowledge of 
the initial properties of the initial black holes is of great importance in several fields. In 
astrophysics, it provides information on the properties of isolated stellar-mass black holes 
produced at the end of the evolution of a binary system of massive stars. In cosmology, it 
can be used to model the distribution of masses and spins of the supermassive black holes 
produced through the merger of galaxies ll67l l68l and whether massive galaxies should be 
expected to host a massive black hole at their centre (but see also refs. 16911701 ). In addition, 
in gravitational-wave astronomy, the a-priori knowledge of the final spin can help the detection 
of the ringdown ll26l . 

I have reviewed those semi-analytic approaches to model the final state from binary black- 
hole coalescences and that have been derived to bypass the still very expensive numerical- 
relativity calculations and obtain a complete and (reasonably) accurate description of the full 
space of parameters. Although following different routes and approximations, the several 
approaches proposed so far have provided a rather accurate description of the final state from 
binary coalescences. For the first time after decades, these approaches and the numerical- 
relativity calculations behind them, have placed us in the position of knowing the properties 
of the final black hole. However, our understanding of these properties has not progressed 
equally rapidly. It is now possible to compute the second and third-order dependences on 
the initial spins and mass ratios but many of these dependences are still unclear. Because 
knowledge without understanding is often sterile, the works reported here should serve as 
stimulus to other studies where the numerical results represent the "background" to be used 
by perturbative or PN studies to go from knowledge to understanding. The recent works in 
refs. l23l [37l l39l l6l 17 1 1 l72l are encouraging and useful steps in this direction. 
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